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Abstract 

Given a planar triangulation, a 3-orientation is an orientation of tlie in- 
ternal edges so all internal vertices have out-degree three. Each 3-orientation 
gives rise to a unique edge coloring known as a Schnyder wood that has proven 
powerful for various computing and combinatorics applications. We consider 
natural Markov chains for sampling uniformly from the set of 3-orientations. 
First, we study a "triangle-reversing" chain on the space of 3-orientations 
of a fixed triangulation that reverses the orientation of the edges around a 
triangle in each move. It was shown previously that this chain connects the 
state space and wc show that (i) when restricted to planar triangulations 
of maximum degree six, the Markov chain is rapidly mixing, and (ii) there 
exists a triangulation with high degree on which this Markov chain mixes 
slowly. Next, we consider an "edge-flipping" chain on the larger state space 
consisting of 3-orientations of all planar triangulations on a fixed number of 
vertices. It was also shown previously that this chain connects the state space 
and we prove that the chain is always rapidly mixing. The triangle-reversing 
and edge-flipping Markov chains both arise in the context of sampling other 
combinatorial structures, such as Eulerian orientations and triangulations of 
planar pointsets, so our results here may shed light on the mixing rate of 
these related chains as well. 

1 Introduction 

The 3-orientations of a graph have given rise to beautiful combinatorics and compu- 
tational applications. A 3-orientation of a planar triangulation is an orientation of 
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the internal edges of the triangulation such that every internal vertex has out-degree 
three. In this paper we study natural Markov chains for sampling 3-orientations 
in two contexts, when the triangulation is fixed and when we consider the union 
of all planar triangulations on a fixed number of vertices. In the case that the tri- 
angulation is fixed, we consider moves that reverse the orientation of edges around 
a triangle if they form a directed cycle. We show that the chain is rapidly mixing 
(converging in polynomial time to equilibrium) if the maximum degree of the trian- 
gulation is six, but can be slowly mixing (requiring exponential time) if the degrees 
are unbounded. To sample from the set of all 3-orientations of triangulations with 
n vertices we use a simple "edge-flipping" chain and show it is always rapidly mix- 
ing. These chains arise in other contexts such as sampling Eulerian orientations 
and triangulations of fixed planar point sets, so there is additional motivation for 
understanding their convergence rates. 

1.1 3-orientations 

Given an undirected graph G = (V, E) and a function a : V" — > Z+, an a-orientation 
is an orientation of E where each vertex v has outdegree a{v). Several fundamental 
combinatorial structures - spanning trees, bipartite perfect matchings, Eulerian ori- 
entations, etc. - can be seen as special instances of a-orientations of planar graphs. 
We refer the reader to [12 [13 [TS] and references therein for extensive literature 
on the subject. Not surprisingly, counting a-orientations is ^P-complete. Namely, 
consider an undirected Eulerian graph G (with all even degrees) ; the a-orientations 
of G, where a{v) = d{v) /2, correspond precisely to Eulerian orientations of G. The 
latter problem has been shown to be #P-complete by Mihail and Winkler [S^ , and 
more recently Creed ^ showed that it remains #P-complete even when restricted 
to the class of planar graphs. 

The term orientation refers to an a-orientation of a planar triangulation where 
all internal vertices (vertices not bounding the infinite face) have a{v) — 3 and all 
external vertices have a{v) — 0. Each 3-orientation gives rise to a unique edge col- 
oring, known as a Schnyder wood, whose many combinatorial applications include 
graph drawing |26l [7] and poset dimension theory [271 . Several intriguing enumer- 
ation problems remain open, such as the complexity of enumerating 3-orientations 
of a planar triangulation (see e.g., [15..) We study the problem of sampling 3- 
orientations of a fixed (planar) triangulation, as well as sampling 3-orientations of 
all triangulations with n internal vertices. In particular, we analyze the mixing 
times of two natural Markov chains for these problems, which were introduced 
previously but had thus far resisted analysis. 

1.2 Results 

First, we study the problem of sampling 3-orientations of a fixed triangulation, 
which was stated as an open problem by Felsner and Zickfeld [TS]. Although 
there is no known efficient method for counting exactly, there are polynomial- 
time algorithms for approximately counting and sampling 3-orientations due to a 
bijection with perfect matchings of a particular bipartite graph (see Section 6.2 in 
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[T5]). This bijection allows us to sample 3-orientations in time 0*{'n7) using an 
algorithm due to Bezakova et al. |lj (improving on the celebrated results of Jerrum, 
Sinclair and Vigoda [E]), but this approach is indirect and intricate. 

We consider instead a natural "triangle-reversing" Markov chain, Mtr, that 
reverses the orientation of a directed triangle in each step, thus maintaining the 
outdegree at each vertex. Brehm [B] showed that for any fixed triangulation T, 
Mtr connects the state space \1/(T) of all 3-orientations of T. We also consider 
a related "cycle-reversing" chain, Mcr, that can also reverse directed cycles con- 
taining more than one triangle. The chain Mcr is a non-local version of Mtr 
based on "tower moves" reminiscent of those in [5D]. We show that both of these 
chains are rapidly mixing. Let A/(T) denote the maximum degree of any internal 
vertex of T. We show: 

Theorem 1. If T is a planar triangulation with Ai{T) < 6, then the mixing time 
of Mcr on the state space "^{T) satisfies 

r(e) = 0(n5 In £-1). 

We can use a standard comparison argument together with Theorem [l] to infer a 
bound on the mixing time of the triangle-reversing chain Mtr- Thus we prove: 

Theorem 2. If T is a planar triangulation with Aj{T) < 6, then the mixing time 
of Mtr on the state space ^'(T) satisfies 

T(e) = 0(n«lne-i). 

Note that the class of planar triangulations with A/ < 6 is exponentially large in n, 
the number of vertices. An interesting special case is when the fixed triangulation 
is a finite region A of the triangular lattice, since sampling 3-orientations on A 
corresponds to sampling Eulerian orientations. Creed [H] independently solved the 
sampling problem in this special case using an approach similar to ours based 
on towers; he shows that for certain subsets of the triangular lattice the tower 
chain can be shown to mix in time 0(n'*). In addition, it was previously shown 
that similar cycle-reversing chains are rapidly mixing in the context of sampling 
Eulerian orientations on the Cartesian lattice ^20J and the 8- vertex model pjj . Our 
analysis here bounding the mixing rate of Mcr. in the general setting of arbitrary 
planar graphs with maximum degree 6 requires additional combinatorial insights 
because we no longer have the regular lattice structure. In particular, we make use 
of a combinatorial structure outlined by Brehm [6]. In fact, this structure allows 
us to extend our analysis to certain non-4-connected triangulations that can have 
vertices of degree greater than six. 

Next, we prove that when the maximum degree is unbounded, the chain Mtr 
may require exponential time. Specifically, we prove: 

Theorem 3. For any (large) n, there exists a triangulation T of size n for which 
the mixing time of Mtr. on the state space ^(T) satisfies 

T(e) = ri(2"/4). 
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Based on the construction we give here, Felsner and Heldt [T3| recently constructed 
another, somewhat simpler, family of graphs for which the mixing rate of Mtb. and 
McR is exponentially large. However, we note that their family also has maximum 
degree that grows with n. It would be interesting to know if there is such a family 
with bounded degree. 

The second problem we study is sampling from 5'„, the set of all 3-orientations 
arising from all possible triangulations on n internal vertices. The set is known 
to be in 1-1 correspondence with all pairs of non-crossing Dyck paths, and as 
such has size C„+2C„ — C^+i, where C„ is the nth Catalan number. Since exact 
enumeration is possible, we can sample using the reduction to counting; this was 
explicitly worked out by Bonichon and Mosbah [S] . We consider a natural Markov 
chain approach for sampling that in each step selects a quadrangle at random, 
removes the interior edge, and replaces it with the other diagonal in such a way as 
to restore the out degree at each vertex. Bonichon, Le Saec and Mosbah [4] showed 
the the chain Mef connects the state space 'i'n and we present the first bounds 
showing that the chain is rapidly mixing. Although the exact counting approach 
already yields a fast approach to sampling, the chain Mef compelling because 
it arises in other contexts where we do not have methods to count exactly. For 
example, it has been proposed as a method for sampling triangulations of a fixed 
planar point set, a problem that has been open for over twenty years. In addition, 
there is additional interest in the mixing rate of this chain precisely because the 
number is related to the Catalan numbers; there has been extensive work trying 
to bound mixing rates of natural Markov chains for various families of Catalan 
structures (see, e.g., [H]). Specifically, we prove: 

Theorem 4. The mixing time of A4ef on the state space '^n satisfies 

T(e) = 0(ni°log(n/e)). 

1.3 Techniques 

The primary challenge behind the proofs of these results is extracting the right com- 
binatorial insights to understand the dynamics in the context of Schnyder woods 
and 3-orientations. Fortunately, there is a long history examining the rich structure 
of Schnyder woods. We were able to extend these results in several new ways, thus 
allowing us to bound the mixing rates of these chains. The proof of Theorem [T] 
for McR. involves a complex coupling argument that is straightforward if T is the 
triangular lattice, but requires more work to generalize to all triangulations with 
Aj < 6. To prove Theorem [3] we produce an intricate triangulation T which is 
carefully constructed to reveal an exponentially small cut in the state space ^'(r). 
Although our choice of T may seem complicated, it was carefully architected using 
properties of 3-orientations to show that the Markov chain may be slow. Our proof 
of Theorem [4] involves a detailed application of the comparison method to bound 
the mixing time oi A4ef by relating it to a local Markov chain on Dyck paths, 
AioK, whose mixing time is known (see [5D1 [H]). The key obstacle here is de- 
composing moves oi Mdk into moves oi Mef while avoiding congestion. This is 
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especially challenging because although M dk is local in the setting of Dyck paths, 
in the context of 3-orientations it can make global changes to a 3-orientation in a 
single step. 



2 Preliminaries 

We begin with background on 3-oricntations, Schnyder woods, and Markov chains. 
The external vertices and edges are those on the infinite face of a planar trian- 
gulation. All other vertices and edges are called internal. Let Tn be the set of 
planar triangulations with n internal vertices. Given a triangulation T e 7^, let 
^'(T) be the set of 3-orientations of T, where a 3-orientation is an orientation of 
the internal edges of T such that every internal vertex has out-degree 3 and ev- 
ery external vertex has out-degree 0. Let — ^tgt„.'^{T) denote the set of all 
3-orientations over all triangulations on n internal vertices. There are several inter- 
esting bijections between 3-orientations of triangulations and other combinatorial 



objects; for example, in Section 4.1 we present a bijection between and pairs 
of Dyck paths. Fraysseix and Ossona de Mendez defined a bijection between vlf(T) 
and the Schnyder woods of T [16] . A Schnyder wood (see Figure |9| is a 3-coloring 
and orientation of the edges of T such that for every internal vertex v, 

• V has out-degree exactly 1 in each of the 3 colors blue, red and green, and 

• the clockwise order of the edges incident to v is: outgoing green, incoming 
blue, outgoing red, incoming green, outgoing blue and incoming red (see 



Figure la) 




Figure 1: (a) The vertex condition, (b) An example of a 3-orientation of a trian- 
gulation on 9 internal vertices. 



In our figures, we differentiate the colors of edges in the Schnyder woods by 
dashed lines (green), dotted lines (red), and solid lines (blue). See Figure lb for an 
example of a Schnyder wood. The orientation of the edges of a Schnyder wood is a 
3-orientation and each of the colors forms a directed tree which spans the internal 
vertices and is rooted at one of the external vertices. See [IS] for a proof of the 
bijection. 
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It is often convenient to consider the Schnyder wood associated with a given 
3-orientation, because the vertex condition allows us to infer some of the global 
information about the 3-orientation by looking locally at each vertex. Throughout 
the paper, when we refer to the colors of the edges of a 3-orientation, we mean the 
colors of the Schnyder wood associated with that 3-orientation. We will use the 
additional information provided by the bijection with Schnyder woods extensively 
throughout the proofs. 

Next, we present some background on Markov chains. The time a Markov 
chain takes to converge to its stationary distribution tt is measured in terms of the 
distance between tt and V*, the distribution at time t. The total variation distance 
at time t is 

\\r\n\\t, = "fiilT. \V\x,y)-7r{y)l 
ye* 

where V^{x,y) is the t-step transition probability and is the state space. For all 
e > 0, the mixing time r of is defined as 

r(e) = min{< : H'P*', 7r|lf„ < e,Vi' > t}. 

We say that a Markov chain is rapidly mixing if the mixing time is bounded above 
by a polynomial in n, where in this case, n is the number of internal vertices of the 
triangulations. 



3 Sampling 3-orientations of a fixed triangulation 

First we consider a Markov chain for sampling the 3-orientations of a given trian- 
gulation. Let T be a planar triangulation with n internal vertices. Consider the 
following natural local Markov chain M.tr on the set of all 3-orientations of T. 
Select a directed 3-cycle at random and reverse its orientation. We will see that 
AAtr samples from the uniform distribution, but its efficiency will depend on T . 



In Section 3.2 we show that if the maximum degree of internal vertices is at most 6, 



M.TR is rapidly mixing. In contrast, in Section 3.3 we demonstrate a triangulation 
T (with unbounded degree) for which A4tr takes exponential time to sample from 
the state space ^'(T). Define Mtr. as follows (see Figure [2]). 

The Markov chain Mtr 

Starting at any (ToG5'(r), iterate the following: 

- Choose a triangle t in ai u.a.r.[^ 

- If i is a directed cycle, with prob . ^ reverse t to obtain cr^+i. 

- Otherwise, (Ti+i = Ui . 

Brehm proved M.tr connects the state space ^'(T) Since all valid moves have 
the same transition probabilities, A4tr converges to the uniform distribution over 
the state space \['(T). 



^The abbreviation u.a.r. stands for uniformly at random. 
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(a) (b) 
Figure 2: A move of M.tr- The triangle Aa6c reverses to Aacb. 

3.1 Background on 3-orientations of Planar Triangulations 

In this section we will provide an overview of several results on 3-orientations of 
planar triangulations and the Markov chain Mtr which we will use in Section 3.2 to 
show that Mtr converges quickly when the maximum degree of the triangulation 
is at most 6. Brehm [6] provides a detailed analysis of the robust structure of 
3-orientations of planar graphs. In particular, he constructs a framework which 
shows that the set of 3-orientations form a distributive lattice and that for a planar 
triangular graph G, any two 3-orientations of G are connected by a series of moves 
of the Markov chain Mtr- As part of this effort, Brehm examines a potential 
function on the faces of the graph, which we will see is useful to upper-bound 
the number of 3-orientations for a given triangulation and the maximum distance 
between any two 3-orientations. 

To show connectivity of Mtr, Brehm considers first the case of 4-connected 
planar triangulations. In this case, every triangle of G is a face and he shows that 
it is possible to get between any two 3-orientations by a sequence of rotations of 
directed facial triangles. Suppose now that a planar triangulation G has exactly 
one triangle t which does not bound a face. For any non-facial cycle, Brehm shows 
that in any 3-orientation of G, the edges in the region bounded by that cycle that 
are incident to some vertex v of the cycle must be directed towards v. This implies 
that no face / contained within t that shares an edge with t can be bounded by a 
directed triangle; hence such faces can never be rotated. In fact, this implies (see [6] 
for details) that G can be regarded as the cross product of the triangulation G\t, 
the restriction of G to the vertices on the boundary of t and the vertices contained 
within the region bounded by t and the triangulation G\t obtained by removing 
all vertices and edges contained within the region bounded by t. Thus by allowing 
Mtr. to rotate arbitrary directed triangles, this amounts to extending Mtr to the 
triangulation G\t which is now 4-connected. The same arguments will hold when 
G has many non-facial triangles. Thus he obtains the following theorem. 

Theorem 5 (Brehm). For any planar triangulation G, Mtr connects the set of 
all 3-orientations ofG. 

In our setting, we use the fact that G\t is independent of G \ i to show that the 
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mixing time of Mtr is the maximum of the mixing times of each 4-connected piece 
of G, subject to the delay which results from the fact that Mtr only attempts to 
update one 4-connected piece at a time. 

Brehm defines a potential X of a 4-connected planar triangulation as follows. 

Definition 1. A potential X of a 4-connected planar triangulation G is a mapping 
f Xf from the interior faces to the natural numbers such that 

• Xf ~ if the boundary of f contains an exterior edge 

• \xf — Xg\ < 1 holds for any two adjacent faces f,g. 

The value of a potential X is defined by Xf. 

It turns out that there is a bijection between 3-orientations of G and a subset of 
the potentials of G, called induced potentials, and that each move of Mtr changes 
the potential of a face by ±1 (see [5] for details). Every triangulation G has 
2n+ 1 faces (not counting the infinite face). It is easy to see that for any potential 
of G, the maximum value for any face is at most [ ^"^^ J since each face can only 
differ from its neighbors by at most 1 and faces adjacent to the boundary have 
value 0. Thus, the maximum number of steps of M.tr required to get between 
two 3-orientations of G is bounded by the maximum value of any potential of G. 
This implies Lemma [6]ja). Moreover, the number of 3-orientations of a graph G is 
bounded by the number of induced potentials of G. Since each face in a potential 
is within 1 from each of its adjacent faces, the number of induced potentials is at 
most 3^"+!. 

Lemma 6. Let G be a ^-connected planar triangulation. 

(a) The maximum distance between two 3-orientations ofG is at most {2n + 1)^/2. 

(b) The number of 3-orientations of G is at most 3^"^^. 

3.2 Mtr mixes rapidly for A/(r) < 6 

In this section we prove Mtr is rapidly mixing on the state space ^'(r), if T is 
a planar triangulation with A/(r) < 6. First, we introduce an auxilliary chain 
McRi which we will then we use to derive a bound on the mixing time of Mtr- 
The Markov chain Mcr involves towers of moves of Mtr, based on the nonlocal 
chain introduced in [20] . Notice that if a face / cannot move then two of its edges 
have the same orientation and the other edge does not (see, e.g., the face /i in 
Figure [3^). We call this edge the disagreeing edge of /. Define a tower of length 
/c to be a path of faces fi, f2, ■ ■ ■ , fk such that the following three conditions are 
met: fk is the only face which is bounded by a directed cycle (i.e. it has a move); 
for every 1 < i < fc, the disagreeing edge of /; is also incident to /i+i; and every 
vertex v is incident to at most three consecutive faces in the path (see Figure [3]) . 
The idea is that once the edges of fk are reversed, then the edges of fk-i can be 
reversed, and so on. We call /i the beginning of the tower, and fk the end. Notice 
that every face is the beginning of at most one tower (it may be a tower of length 



8 



1). The effect of making these moves is to reverse the edges of the directed cycle 
surrounding the tower (although the colors on the internal edges also change). 



The Tower Markov chain J^cR 



Starting at Einy ctq , iterate the following: 

• Choose a (finite) face / u.a.r. 

• If / is the beginning of a tower of length k, then 

k>2 



w.p. 

• Otherwise, Ui+i = tii . 



k = 1 



reverse this tower to obtain ai+i. 



The moves of A4tr are a subset of the moves of A4cb., so AicR is connected 
as well. 




Figure 3: A tower of length 6. 



We first consider the case that T is 4-connected. Notice if T is 4-connected, 
every 3-cycle is facial, so M.tr selects a face and rotates the edges around that 
face if possible. The bulk of the work to prove Theorems [l] and |2] is to show that 
M.CR is rapidly mixing when T is 4-connected. 

Theorem 7. Let T he a ^-connected planar triangulation with A/(T) < 6. Then 
the mixing time of Ai cR, on the state space ^!{T) satisfies 

r(e) = 0(n5 In £-1). 

First we prove Theorem [7] and then we apply the comparison theorem to prove 
that Mtr is also rapidly mixing and extend the result to non-4-connected trian- 
gulations using a result of Brehm [B] , proving Theorems [l] and [2) The main tool we 
will use to prove Theorem [7] is the following path coupling theorem, due to Dyer 
and Greenhih [T^. 

Theorem 8 (Dyer and Greenhill). Let (j) he an integer-valued metric on ^'x^ which 
takes values m {0, B}. Let U he a suhset of'^x'^ such that for all {xt,yt) G ^'x^E' 
there exists a path xt — zq, zi, Zr — yt between xt and yt such that {zi, Zi+i) G U 
for < i < r and J2l=o ^i+i) = 'Pi^t, Vt)- Let Ai be a Markov chain on state 
space 5" with transition matrix P. Consider any random function / : 5" — > ^ such 
that P[f{x) — y] = P{x,y) for all x,y ^ , and define a coupling of the Markov 
chain by {xt,yt) {xt+i,yt+i) = {f{xt),f{yt))- If E[(j){xt+i,yt+i)] < (j){xt,yt), 
for all Xt,yt & U, let a > satisfy FT[4>{xt+i,yt+i) 4>{xt,yt)] > ce for all t such 
that Xt yt- Then the mixing time of A4 on the state space ^> satisfies 

eB^' 



T(e) < 2 



[Ine-^l. 
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Proof of Theorem^ Let T be a 4-connected planar triangulation with A/(T) < 6. 
Define the distance d between any two 3-orientations in the state space 'if{T) to be 
the minimum number of steps of Mtr from one to the other. Assume cr, r G 'J'(T) 
and r is obtained from a by reversing a facial triangle /. We use the trivial 
coupling, which chooses the same face for a and r at every step. Suppose without 
loss of generality that the edges of / are clockwise in cr. There are two obvious 
moves that decrease the distance, namely when the AicR selects the face / and 
chooses to direct the cycle clockwise or counterclockwise, each of which happens 
with probability l/(2(2n + 1)). Moreover, any move of Mcr that does not involve 
an edge of / occurs with the same probability in a and r, and hence is neutral. 

We call a tower bad if it contains a neighbor /' of / that is not the end of 
the tower. In this case, we say this bad tower is associated with f . On the other 
hand, a tower is good if it ends in /, or if it ends in a face /' adjacent to / and 
contains no other faces adjacent to /. We will show that the good towers in a have 
corresponding good towers in r, while the bad towers in a fail in r, and therefore 
increase the distance. Any tower that is neither good nor bad does not contain an 
edge of /, so it is neutral with respect to the distance. 

Suppose fc > 1, (/i, /2, . . . , fk) is a good tower in cr, and fk is adjacent to /. 
We claim that (/i, /2, . . . , /fc, /) is a good tower in r. It is clear that in r, / is the 
only one of these faces that is bounded by a cycle, and that upon rotating /, the 
tower (/i, /2, • • ■ , /fe) is possible. We must check two things: that (/i, /2, ■ • ■ , /fc, /) 
is a path of faces (i.e. does not contain any cycle of faces), and that every vertex 
is incident to at most three consecutive faces. The first condition is clear, since fk 
is the only neighbor of / in {/i, /2, . . . , fk}, and (/i, /2, ■ • ■ , /fc) is a path of faces. 
Suppose the second condition does not hold. Then there is a vertex v incident 
to f,fkjfk-i, and /fc_2. The edges between faces /fe_2 and fk-i and between 
fk-i and fk are either both incoming to v or both outgoing from v (see Figure |4]). 
Moreover, since the edge between fk-2 and fk-i is the disagreeing edge of /fe-2, 
the two edges of /fe_2 incident to v are either both incoming to v or both outgoing 
from V (similarly the two edges of fk incident to v are either both incoming to v or 
both outgoing) . Hence there are four edges incident to v which are all incoming or 
all outgoing; a contradiction since a vertex of degree at most 6 with exactly three 
outgoing edges can have at most three incoming edges as well. Therefore if a good 
tower of length fc > 1 begins on a face /i and ends on a neighbor fk of / in ct then 
there is a corresponding tower of length fc + 1 that begins on /i and ends on / in 

T. 

Thus we have shown that if (/i, /2, • ■ • , fk) is a good tower in cr, then (/i, /2, . . . 
...,fk,f) is a good tower in r. On the other hand, it should be clear that if 
ifi, f2, ■ ■ ■ , fk, f) is a good tower of length fc + 1 > 2 that ends on / in cr then 
ifi, f2, ■ ■ ■ , fk) is a good tower of length fc in r. In either case, if fc > 2 then the 
expected change in distance given the choice of these towers is 

1 ( 1 i/Yl 1 

2n + l V 6(fc + l) \6k 6(fc + l) 
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(a) (b) 

Figure 4: A tower (/i, /2, . . . , fk) in (J for which (/i, /2, • • • , A, /) is not a tower in 
T. Notice deg(w) > 7. 



If fc = 2 then the expected change in distance is 

1/1 /I 1\\ 1 



2n + lV 12 V2 12// 3(2n + l)' 

We point out that if cr and r have good towers using a neighbor /' of / then 
no bad tower in cr or r is associated with /'; that is, if there exists a bad tower 
containing /' then /' is the end of the tower. Suppose without loss of generahty 
that the good tower is longer in r than in a. Then the edge between / and /' 
is the disagreeing edge of /' in a so the only way to tower is towards /, so /' is 
not in a bad tower in a. On the other hand, in r, /' is bounded by a cycle, so 
it must be the end of any tower containing it. We claim that cr and r can have 
at most two bad towers associated with a given face /' adjacent to /. It is clear 
that a (resp., r), has at most one bad tower that begins in /', which is defined by 
the disagreeing edge of /' in a. However, a may have a bad tower that uses /' 
but does not begin in /'. Let (/i, /2, • • ■ , /fe) be such a tower. We will show that 
/' = f2', suppose not, so that /' = /i, where j > 3. Then as above, there is a vertex 
V that is incident to /, fi, fi-i, and fi-2, and the same proof will show that v must 
have either in-degree at least 4 or out-degree at least 4, which is a contradiction. 
Therefore bad towers associated with /' must either begin in /' or in a neighbor 
/i of /'. Moreover, if there is a bad tower in cr (r) beginning in fi then in both 
cr and r, the edge between faces /' and /i is /I's disagreeing edge, which means 
that r (resp. cr) cannot have a bad tower beginning in /'. Therefore there are at 
most two bad towers in cr or r associated with /'. The expected change in distance 
given that a bad tower of length fc > 2 is chosen is 2fc/(6A:(2n + 1)) = l/(3(2n-|-l)). 
Therefore 

1 / . /1\ . /I' 



^[^^]^2^l-H2j+H3jj=°- 

We will see that for any {(Jt,Tt) (not necessarily at distance 1), Pr (c?(crt, Tt) ^ 
d{at+i,Tt+i)) > l/(4n + 2). This follows from the connectivity proof of Brehm [5]. 
For any such pair, there exists a path of transitions of Mtr from at to Tt of 
length d{at,Tt), and the first of these transitions occurs with probability at least 
l/(4n-|-2), since each face is selected with probability at least l/(2n-|- 1), and the 
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move succeeds with probability 1/2. Clearly this move decreases the distance (the 
number of steps between the configurations) by 1. By Theorem [s] and the bound on 
the distance between any two 3-orientations given in Lemma[6^, we see the mixing 
time of M.CR over 4-connected triangulations satisfies 

r(e) < \e{2n + 1)^21 [Ine-^l = 0{n^ Ine'^). 

□ 

Next, we prove Theorems [T] and [2] We will use Theorem [7] and the comparison 
method, Theorem [T2| introduced in Section [4?2| to bound the mixing time oi Mtr 
in terms of the mixing time of McR in the case of 4-conncctcd triangulations. We 
will also extend the analysis to all planar triangulations by showing that Mtr. 
operates as a product of independent Markov chains, each acting on a 4-connected 
planar triangulation. Thus, we will need one final detail, which is the following 
straightforward theorem, proved in [2^ . 

Theorem 9. Suppose the Markov chain M is a product of M independent Markov 
chains Mi, A42, ■ ■ ■ tM.m, where M updates Mi with probability pi, where pi = 
1. If Ti{e) is the mixing time for Mi on the state space 4* and Ti(e) > line for 
each i, then the mixing time of M on the state space ^ satisfies 

rfe) < max — r, ( ^ . 

^ ' - z=l,2,...,M p, \2MJ 



Proof of Theorems^ and^ 

First, we compare Mcr with Mtr using Theorem [T2| to derive a bound on 
the mixing time of Mtr in the case of 4-connected planar triangulations. To do 
so we need to bound the constant A given in that theorem. For each edge {a, r) 
in Mtr that takes a counterclockwise cycle / and makes it clockwise, let T{a, r) 
denote the set of edges (x, y) of Mcr such that the tower /i, /2, • ■ • , /fe that takes 
X to y contains the face /. Since a tower is uniquely determined by its direction 
and its length, |r(cr,T)| < 3n^. Define A :— max(g..r) (|r(CT, r)|) < 3n^. Since 
each 3-orientation has the same stationary probability. Lemma [6]d implie s th at the 
minimum weight of any state is tt* > 3-(2n+i) Therefore, by Theorems 12 and[7j 
the mixing time of Mtr over 4-connected planar triangulations satisfies 

log( — ) 

Tie) < 4 At' 

log(l/(26)) 



log(l/(2e)) 



O I ,f ' ■ nMne-i 



/ n-loge 7 _i 
Vlog(l/(2e)) 
= 0(n^ Ine^i). 

Finally, we can extend this to non-4-connected planar triangulations, where 
Mtr may select non-facial triangles. Brehm [B] proves that if T has a non-facial 
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triangle C, the edges on its interior that are incident to C must point towards C. 
This implies that for all a G \I'(r) every face on the interior of C that contains an 
edge of C is not bounded by a directed cycle, so cannot rotate, regardless of the 
orientation of C . Thus Mtr acts completely independently on the interior of C 
from the exterior of C, and its mixing time is independent of the mixing time of 
the exterior as well. 

Let T be a planar triangulation and assume that Ci , C2 , . . . , Cfe are all the non- 
facial triangles of T. Let Ti be the triangulation consisting of all faces contained 
within Ci and not within any other non- facial triangle contained within Ci. Let Ti 
be the mixing time of Mtr on Ti,Fi be the number of faces within Ci, and rij be 
the number of internal vertices to Ti . Then n = ^ ^ rii and the number of faces in 
T is Fi. Therefore, by Theorem |9j the mixing time of A4tr on T will be 



r(e) = max 



e \ 2n + l 
2k) F, 

2k\ 2n + 



O ( nf In 



This is maximized when ni = n, k = 1, and i^i = 2n + 1, so 

T(e) = 0(n«ln(e-i)). 

This proves Theorem [2] In fact, this shows that A4tr is rapidly mixing on the 
state space \['(r) for any planar triangulation T whose 4-connected triangulations 
Ti, T2, . . . , Tfe each have maximum degree (of any internal vertex) 6. Moreover, the 
same proof can be applied to the bound on A4cR for 4-connected triangulations to 
extend to the case of non-4-connected triangulations, proving Theorem [T] 

□ 



3.3 M.TR can mix slowly 

We now exhibit a triangulation on which Mtr takes exponential time to converge. 
A key tool is conductance, which for an ergodic Markov chain with distribution tt, 
is $ = min sc* X^siss "'('^i^^^'^i' '^2)/7r(S'). The following theorem relates 

7r(S)<l/2 ' ___ 

the conductance and mixing time (see [13 )• 

Theorem 10. For any Markov chain A4 with conductance the mixing time of 
A4 on state space 5* satisfies 

r(e) > (4$)-i - 1/2. 

Proof of Theorem For the generalized triangulation G given in Figure [5] with 
n = 4t — 2 internal vertices, Mtr takes exponential time to converge. Specifically, 
we show that although there is an exponential number of 3-orientations where 
the edge (uo,Wt+i) is colored blue or red, all paths between these 3-orientations 
with {vo,vt+i) colored differently must include a 3-orientation where {vo,vt+i) is 
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Figure 5: A triangulation for which AAtr mixes slowly. 



colored green. Moreover, we show there is only a single 3-orientation that satisfies 
this property (namely, the one pictured in Figure |5|, which creates a bottleneck in 
the state space. 

Let D be the set of 3-orientations of G with (wq, Vt+i) colored red or green and 
Z?, be the set of 3-orientations with [vq^vi+i) colored blue. In order to show that 
both D and D are exponentially large we produce a triangulation in each set which 
contains roughly t directed triangles which do not share any edges and reversing 
these triangles does not change the colors of the edges adjacent to vq. Hence each 
of the 2* choices of the orientations of these triangles gives a distinct 3-orientation 
with edge {vQ,Vt+i) colored appropriated. 

First consider the 3-orientation in Figure [6] which is in D. Notice that trian- 
gles Ti,T2, ....Tt-2 do not share any edges and reversing these triangles does not 
change the color of any edges adjacent to wq- Each of these triangles has 2 possible 
orientations and each of these 2* choices of the orientations of the triangles gives a 
distinct 3-orientation with edge (wojWt+i) colored red implying that 

\D\ > 2* = 2("-6)/^ (1) 

Next, consider the coloring in Figure [7] Notice that triangles S'l, 5*2, ....S"*-! do 
not share any edges and reversing these triangles does not change the color of any 
edges adjacent to vq. Each of these triangles has 2 possible orientations and each of 
these 2*~^ choices of the orientations of the triangles gives a distinct 3-orientation 
with edge (fo,Wj+i) colored blue implying that 

\D\ > 2*-i = 2("-2)/'^. (2) 
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Figure 6: There is an exponential number of 3-orientations with edge (vQ,Vt+i) 
colored red corresponding to the different orientations of triangles Ti,T2,. . . , Tt-2- 

Next, we show that there is only one 3-orientation of G with [vo, wt+i) colored 
green, corresponding to Figure [Sj By the Vertex Condition, if edge {vo^vt+i) is 
green then edges (wq, wi), (wo, 1^2), (fo, I't) rnust all be directed toward vq and 
colored red; this is because the edge {vq, suue) is blue and directed towards Sbiue in 
every 3-orientation of G. Similarly, edges (vq, wt+2), (vq, vt+3), . . . , (vq, V2t+i) must 
all be blue and directed toward vq. Since Wj+i has degree 4 and has the incoming 
green edge {vo,Vt+i), the other edges incident to Vt+i are determined; {vt,Vt+i) 
is blue, (vt+i,V3t+i) is green and (wf+i, 114+2) is red all directed away from vt+i- 
Knowing the colors and direction of the edges incident to vq and vt+i forces the 
color and direction of the edges incident to Vt+2- Similarly vt+i and vt+2 force v^t+i 
and Wo I 'Vt+2 and v^t+i force Vt+3 and so on until the color and direction of all edges 
incident to vertices Wt+i, Vt+2 • ■ • , ^24+1 and V3t+i,V3t+2, ■ • ■ , ^^44-2 are forced. Next 
consider Vt] we know the edges (w*, t't+i), (■^t, ^^3t+i), ('^t, i'3t+2), • ■ • , f4t-2) are 
all blue and directed toward Vt which implies (w*, wat) and Vt-i) must be directed 
outward and green and blue respectively. Now consider v^t ; we have already shown 
that all edges incident to v^t except for {vstjVt-i), {v3t,V3t-i), {v3t,Vsi) have a 
forced color and are directed inward. Thus, these 3 edges must all be directed 
outwards with colors blue, green and red respectively. Similarly, knowing the colors 
and directions of all edges incident to v^t , Vt and vq forces the colors and directions 
of edges incident to vt-i and V3t,vt-i forces v^t-i and so on for the remaining 
vertices. Since all the edge colors and orientations are fixed, this implies there is 
a unique 3-orientation with (z)o,Wt+i) colored green. To go from a configuration 
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Figure 7: There is an exponential number of 3-orientations with edge {vQ,Vt+i) 
colored blue corresponding to the different orientations of triangles Si, S2, ■ ■ ■ , St-2- 



where edge (uq, vt+i) has color red (blue) to blue (resp., red) one must go through a 
coloring where the edge is green. This is because the only choices for edge (vq, Vt+i) 
are red directed toward Wg, blue directed toward Vq, and green directed away, and 
any move that changes the color must also change the direction. 

Finally, given the bound on |D| and \D\, we derive a bound on the mixing time 
of Mtr- Let g be the single 3-orientation which has edge {vo,Vt+i) colored green. 
If 7r(D) < 1/2 then combining the definition of conductance with the bound on |D| 
yields 



< 



l^)Ed^eD,d2eD'^idi)r{di,d2) 



i/z 



t(D} Z-^d-ieD KH^"-^! — 7r(D) — 2("-6)/*/Z ~ 2("-6)/'i ' 

If ^{D) > 1/2 then ^{D) < 1/2 and so by detailed balance and the bound on \D\ 
•^Mtb < ^)Ed,eD,d2eD^idi)'P{di,d2) = ^^Ed,eD^idi)'Pidi, 9) 



In both cases, 

'^Mtr < 2^("-6)/4^ Applying Theorem proves that the mixing 
time of A4tr satisfies 



-(e) > 2("-i4)/4) _ I 



□ 
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Remark 1. Combining this result with the comparison argument in Section 3.2 
shows that AicR can also take exponential time to converge. 



4 Sampling the 3-orientations of triangulations on 
n internal vertices 

In this section, we define the local Markov chain Mef for sampling uniformly 
from and show that Mef is always rapidly mixing. Our argument relies on 
a bijection with pairs of Dyck paths to relate the mixing time of a local Markov 
chain on Dyck paths to that of Mef- Define Mef as follows (see Figure [8}d). 

The Markov chain Mef 

Starting at any ctq € , iterate the following: 

- Choose facial triangles Ti and T2 with shared edge xtj u.a.r. 

- Choose an edge from Ti U T2 u.a.r., if one exists. With 
prob. 1/2 replace the path by the path {x^,ziZ;} where w 
is the remaining vertex of Ti U T2 . 

Otherwise, fii+i =cr,;. 



• ^ • • k- 

X \ y X •. 



Figure 8: A red/green swap. 



If the edge with color a is replaced by the edge xh with color Cj , we call this a 
Cj/ci swap. Bonichon, Le Saec and Mosbah showed in [4] that Mef connects the 
state space 5'„. Since all valid moves have the same transition probabilities, this 
implies that Mef converges to the uniform distribution over state space 5'„. 



4.1 The bijection between \l/„ and pairs of Dyck paths 

The key to bounding the mixing time of A^bf is a bijection between and pairs 
of nonoverlapping Dyck paths of length 2n, introduced by Bonichon [3J. Dyck 
paths can be thought of as strings a\a2 - ■ ■ a^n containing an equal number of I's 
and —I's, where for any 1 < fc < 2n, X^iLi ^ 0- Recall that a 3-orientation 
of a triangulation can be viewed as a union of three trees, one in each color. In 
the bijection, the bottom Dyck path corresponds to the blue tree, and the top 
Dyck path indicates the degree of each vertex in the red tree. The green tree 
is determined uniquely by the blue and red trees. More specifically, given a £ 
\E'„, to determine the bottom Dyck path, start at the root of the blue tree and 
trace along the border of the tree in a clockwise direction until you end at the 
root. The first time you encounter a vertex, insert a 1 in the Dyck path, the 
second time you encounter the vertex insert a —1. Let vi,V2, . . . ,Vn be the order 
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of the vertices as they are encountered by performing this DFS traversal of the 
blue tree in a clockwise direction and define L to be the resulting linear order 
on the vertices. Let di be number of incoming red edges incident to Vi. Let 
r be the number of incoming red edges incident to Sred- The top Dyck path 
is as follows l(-l)''n(-l)''n(-l)''^ . . . The structure of the 3- 

orientation guarantees that the top path will never cross below the bottom path. 
See Figure [9) and '3' for details. 




'green 



Figure 9: The bijection between 3-orientations and Dyck paths. 

We will deduce that A4ef is rapidly mixing by comparing it to Mdk, an 
efficient Markov chain on (pairs of) Dyck paths which was introduced by Luby, 
Randall and Sinclair [5D]. The algorithm proceeds as follows. At each step select 
a point on one of the two Dyck paths uniformly at random. If the point is a local 
maximum (or minimum) then push it down (or up) with probability 1/2 as shown 
in Figure [Tor a-b). If this move is blocked by a local maximum (or minimum) in the 



bottom (or top) Dyck path as shown in Figure 10 ; then push both Dyck paths down 
(or up) with probability 1/2 as shown in Figure [To| (c-d). The following theorem 
due to Wilson bounds the mixing time of M dk ■ 



(a) (b) (c) (d) 

Figure 10: Two moves of the Markov chain Mdk- 



Theorem 11 (Wilson). The mixing time oj on state space satisfies 

r(e) = e{n^\og{n/e)). 

Using the above bijection, the Markov chain A4 dk on Dyck paths can be translated 
into a Markov chain on 3-orientations of triangulations, but its moves are quite 
unnatural in that setting. We obtain a bound on the mixing time oi Mef using 
Theorem [11] together with a careful comparison argument. 
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4.2 M.EF mixes rapidly 

Next we show that A^^i? is an efficient algorithm for sampfing from the state space 
\l/„ by comparing the moves of Mef to the moves of Mdk- The comparison 
theorem of Diaconis and Saloff-Coste [5] allows us to compare the mixing times of 
two reversible Markov chains Ai and Ai' on the same state space ^1/. Assume they 
have the same stationary distribution tt and mixing times r and t' respectively. 
Let E(M) = {{X,Y) : M{X,Y) > 0} and E{M') = {{X,Y) : M'iX,Y) > 0} 
denote the transitions of the two Markov chains, viewed as directed graphs. For 
each (X, Y) G E{M'), define a canonical path jxy = {X = Xq, Xi, - ■ ■ ,Xk = Y) 
with {Xi,Xi+i) e M, and let \jxy\ denote the length of the path. Let T{Z, W) = 
{{X,Y) e E{M') : {Z,W) G ^xy} be the set of canonical paths that use the 
transition [Z, W) of M.. Let tt* = minxe* '^{X). Finally, define 

A= max y \-ixY\T^{X)M'{X,Y)/{Tr{Z)M{Z,W)). 

V(ZM) 

We will use the following version of comparison, due to Randall and Tetali |25) . 

Theorem 12 (Randall and Tetali). With the above notation, the mixing time of 
A4 on the state space 4' satisfies 

log( — ) 

Next, we introduce some notation. Let ci be blue, C2 be red, and C3 be green. 
Given a vertex v and i € {1,2,3}, the unique outgoing edge of v with color Ci is 
called v's Ci edge. We also define the first (last) incoming c^-edge of v to be the 
incoming Q-edge of v that is in a facial triangle with u's Ci_i edge (respectively, 
v's Ci_|_i edge, where the subscripts are taken modulo 3). In our canonical paths, 
we will often need to move a Cj edge, say vi, from some neighbor x of w to another 
neighbor y of u across several Ci edges. This is achieved through a sequence of Cj /ci 
swaps as in Figure [TT] 

■\ yy \ jy r. ny r ry 



\ 7,1 



\» / \* , / \> /// 

/ ■■},'/ 

Figure 11: A sequence of red/green swaps. 

Finally, we can prove Theorem [4] 

Proof of Theorem In order to apply Theorem [12] to relate the mixing time of 
Mef with the mixing time of Mdk we need to define, for each transition of 
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A4dk, a canonical path using transitions of J^ef- Then we will show that each 
transition of J^ef is not used by too many canonical paths. We consider four 
cases, depending on whether the move of M.dk takes a peak to a valley and which 
Dyck path it involves. If the move e = (X, Y) affects both paths, we view the move 
as two separate moves {X,Z) and {Z,Y), one on each path, and we concatenate 
the canonical paths as follows: ^x,y — {jx,z,1z.y)- Hence in the following, we 
assume that the transitions of A4 dk affect only the top or the bottom Dyck path. 

A Peak to Valley Move on the Top Dyck Path: Let e = {X,Y) be a 
transition of M.dk that affects the top path only and moves from a peak to a 
valley. Suppose e moves the i"^ 1 (where « > 1) on the top path to the right one 
position (i.e. the Dyck path move swaps the i**^ 1 with a —1 on its right, changing 
a peak to a valley). From the bijection, we know this move does not affect the 
blue tree and corresponds to, in the red tree, increasing the incoming degree of 
Vi by one and decreasing the incoming degree of w^+i by one. If Vi and Wi+i are 
adjacent in the blue tree (there is a blue edge Wi+iWi) then this implies that there 
is a red/green swap involving ViS green edge and Wi+i's first incoming red edge. 
This swap is exactly the peak to valley move, so ^xy = e. Otherwise, we define 
two stages in the canonical path ^xy- 

To assist in defining the canonical paths, let gi be the parent of Vi in the green 
tree. First, we will show that gi is not Sgreen] that is, Vi cannot have its green 
edge point to the green source Sgreen as long as e = {X, Y) is a valid transition 
of J^DK- For the sake of contradiction, suppose ViS green is an edge in X. Let 
P be the path from suue to Vi in the blue tree combined with the edge ViS green- 
This path partitions the vertices of X into two sets according to which side of the 
path they are on. Let S be the set of vertices on the same side of P as Vi+i (not 
including the vertices along P) and let S be the remaining vertices. Let x G V{X) 
and define x's red edge by xrt. We claim that L(x) < L{rx)\ this is a direct result 
of a bijection between Dyck paths and 3-orientations given in [3 (this is a different 
bijection from the one given in Section[2|. Now follow the blue path from x to suue 
until you reach the first vertex z' that has an incoming blue edge from a vertex z 
such that L{z) > L{x). Define z^ = z. We claim that in all 3-orientations with 
this same blue tree, L{rx) > L{zx). To prove this, we need some facts from [3]. 
Specifically, they show that any 3-orientation can be obtained by starting from a 
fixed 3-orientation, called a star realizer, which is unique to each blue tree, and then 
applying a series of transformations. The star realizer satisfies the property that 
L{rx) > L{zx) for any vertex x, and each transformation maintains this property. 

This implies that any red edge in X directed toward vertices in S must be 
directed toward vertices in S for all 3-orientations with the same blue tree, which 
implies that there is no configuration Y which has the same blue tree as X and the 
same incoming red degree sequence except for increasing Vi by one and decreasing 
Vi+i by one. This is a contradiction since we have assumed (X, Y) is a valid move. 
Hence we have shown that gi is not s green- 

Next, let Vj be the parent of gi in the red tree. Notice that L{vi) < L{vi+i) < 
L{vj), since L(vj) > L[gi) and ViS red and green edges prevent Vj from satisfying 
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Figure 12: The vertex Vi's red and green edges prevent vj from satisfying L{vi+i) > 
L{vj) > L{g,). 



Llvi 



> L( 



> L{gi) as shown in Figure 12 In the first stage of the path 



IXY we make the sequence of red/green swaps centered at gi that move the red 
edge giv] to giv\ without affecting any other red edges as shown in Figure 13 d, 
step 1 (see Figure 11 for detail on the sequence of swaps). In the second stage 
we transfer an incoming red edge from w^+i to Vj, completing 7x,y- Recall that 
L{vj) > L{vi+i), so we do this iteratively by moving an incoming red edge 
either to one of x's neighbors in the blue tree that is larger in L or to x's parent in 
the red tree, which is also larger in L, if it is a leaf and has no neighbors as shown 
in Figure [TSJd, step 2 and 3. We claim it is always possible to make one of these 
moves. If X has a neighbor y in the blue tree such that L{y) — L{x) + 1 then there 
must be a green/red swap centered at x and involving x's green edge that moves 
an incoming red edge from x to y as desired. Next, if a; is a leaf with red edge xrt 
then again there is a green/red swap centered at x involving x's green edges that 
moves an incoming red edge from x to as desired. Finally notice that, using this 
canonical path, we never bypass Vj because the original red edge g^Wj blocked any 
blue leaves between w^+i and gi from having red parents higher in L than Vj. 




Figure 13: The canonical path to move a peak down to a valley, top Dyck path. 



Given a transition {Z,W) of A4ef we must bound the number of canonical 
paths 7x1' using this edge. To do so, we analyze the amount of information needed 
in addition to {Z, W) to determine X and Y uniquely. We record the vertex Vi and 
the vertex Vj. If Vi and Vi+i are adjacent we record Ui+i instead of Vj. Notice in 
this case the canonical path only involves red/green and green/red swaps. If we are 
moving a red edge to a higher vertex in L then we are in stage 2 and otherwise we 
are in stage 1. Given this information wc can uniquely recover X and Y . We only 
need to record two vertices, so in this case there are at most canonical paths 
which use any edge {Z, W). 
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A Valley to Peak Move on the Top Dyck Path: Next consider the case 
where e — {X, Y) changes a valley to a peak. Assume that e moves the «th 1 on 
the top Dyck path to the left one position. From the bijection we know that this 
move corresponds to, in the red tree, decreasing the incoming degree of Vi by one 
and increasing the incoming degree of Vi+i by one. If Vi and Vi+i are physically 
adjacent in the blue tree then there is a green/red swap centered at Vi and involving 
Vi's green edge that performs exactly these changes so we will define jxy to be this 
single move. If Vi and Vi+i are not adjacent then we will define two stages in jxy- 
Let Ti be f^'s parent in the red tree. In the first stage we transfer an incoming 
red edge from Vi to using the green/red swap centered at Vi and involving Vi's 
green edge. In the second stage we will repeatedly use the canonical path defined 
for peak to valley moves to transfer an incoming red edge from to Vi+i. 

Again given a transition {Z,W) oi Mef we must upper bound the number of 
canonical paths 'yxY that use this edge in this case. First, we will record a bit to 
determine which stage we are in. For both stages we will record the vertex Vi and 
the vertex r^. If Vi and Vi+i are adjacent we will record vt+i instead of r^. For 
the second stage if we are in the process of making a peak to valley move which 
affects the vertices Vp and Wp+i, we will record the vertex Vp. In addition we will 
record the extra vertex Vj that is needed when we are in the middle of a peak to 
valley move, as explained in the peak to valley case above. Given this information 
we can uniquely recover X and Y. This implies that in this case there are at most 
2?!^ canonical paths which use any edge (Z, W). 





Figure 14: A move of Ai dk that moves the bottom Dyck path from a valley to a 
peak takes the blue tree from (a) to (b). 

A Valley to Peak Move on the Bottom Dyck Path: Let e = (A, Y) be 

a transition of the tower chain McR on the bottom Dyck path, and suppose the 
edge e is a move that takes a valley of the bottom Dyck path and moves it to a 
peak. This affects the blue tree as follows (see Figure [T4| : a's blue edge moves 
from c to b, and all the (blue) children of at: (if any exist) become children of 
be. However, to define the canonical path between these two configurations, it is 
necessary to also know what the top Dyck path looks like, as it determines the red 



(and therefore the green) tree. Figure 15 shows how the red and green tree might 



look. Our path will first update the blue tree from X to match the blue tree of Y, 
and then update the red tree (and therefore the green tree) to match the red tree 
of Y using the steps outlined in the previous two cases. For a vertex u G T let 
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denote the head of v's red edge in X. We will go through 4 distinct stages in the 
canonical path. In stage 1 the blue edge of a moves from c to b. Then in stage 
2, a's red edge moves into position for stage 3, where all incoming blue edges to a 
move down to point to b. Finally, in stage 4 we repair the red tree. 
Stage 1. Given the vertex condition (Figure [la^) for b and the bijection between 
the bottom Dyck path and the blue tree, there is no edge in the angle Zacb and 
thus a, b, and c form a triangle in our triangulation. Vertex b may have some green 
edges coming in between its blue edge be and its red edge bri,. If so, then is 
an edge, as in Figure 15 a,. The first step along the canonical path from X to Y 
is to rotate 6's red edge to point to a (see Figure [Isj^a-b)). This is accomplished 
through a se quence of red/green swaps, one for each green edge coming into b, as 
in Figure Now we are ready to move the blue edge to by swapping with 
b's red edge, as in Figure 15 h-c). 





Figure 15: Canonical path to move a valley up to a peak in the blue tree. 



Stage 2. Consider vertex a. Counterclockwise from a's blue edge, there may be 
some green edges coming into a, followed by the red edge arl. If there are green 
edges, let d be the tail of the last incoming green edge to a; however, if there are no 
green edges, then = d. In this stage of our canonical path, a's red edge moves 
from ra to d (if Tq 7^ d) by a series of red/green swaps. See Figure [TT|c-d). 
Stage 3. Next, we use a's red edge to move all the blue children of alj to point to 



b, one at a time, in a clockwise manner. See Figure 11 'd-e). Now the blue tree is 
completely fixed. 

Stage 4. Finally, we must repair the red tree. Notice that the red edges of a and 
b are the only red edges that we moved; we must move them to their proper place. 
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We must increase the degree of rf, and ra back to match their indegree in X. To do 
this without affecting the blue tree, we first make a's red edge point to ri, (we'h caU 
this Stage 4a) and then make b's red edge point to Va (Stage 4b). These moves 
can each be accomplished by a sequence of red/green swaps without affecting the 
blue tree. 

Given a transition {Z, W) of Mef, we must upper bound the number of canon- 
ical paths jx,Y that use this edge. If {Z, W) is in stage I, we need to remember 
vertices a and r^, and a bit to tell us whether or not we have moved ah yet. Given 
this information, we can recover b and c and can undo all red/green swaps in order 
to get back to X. Given X we can find Y since we know which valley to flip up. 
If {Z,W) is in stage 2, then we only need to record rb, since {Z,W) moves a's 
red edge, so we know a. To get back to the last configuration in stage I, we just 
need to move a's red edge counter-clockwise until it can't make any more red/green 
swaps. Thus we can get back to the last configuration in stage 1, and using a and 
Tb we can recover X. If {Z, W) is in stage 3, we need to record Vh- Each move in 
stage 3 takes a child of m> and moves it to point to b. Hence, we know a. Notice 
that since Aabc was facial in X, all blue edges coming into b in ai before (in 
the counterclockwise direction) were children of a^ in X. Thus, given a, we know 
that we must use a's red edge to move each of these children back up to a. This 
brings us back to the last configuration in stage 2; using a and rj,, we can recover 
X. If {Z,W) is in stage 4a, then we know that the blue tree agrees with Y, and 
the red edges of a and b are the only red edges in a different position in Z than in 
Y. Since {Z, W) moves a's red edge, we know a, and b is the vertex that shares 
a's blue edge. Given b, it is easy to recover ra, since to find r^, just move 6's red 
edge counter-clockwise until it can't make any more red/green swaps. We need to 
record rt and then it is easy to get to Y. If {Z, W) is in stage 4b, then we know 
that the blue tree agrees with Y and 5's red edge is the only red edge that is in 
a different position in Z than in Y. Since {Z, W) moves &'s red edge, we know b. 
As described in stage 4a, it is easy to then find r^; move 6's red edge all the way 
there to get to Y. In each of the four stages we need to record a maximum of 2 
vertices and a single bit. This implies that in this case there are 0{n^) canonical 
paths which use any edge {Z, W). 



A Peak to Valley Move on the Bottom Dyck Path: Next suppose e — 
{X, Y) is a move that takes a peak of the bottom Dyck path and moves it to a 
valley. This moves the blue tree in Figure [Mjj to Figure 14 1. In general, a's blue 
edge moves from b to c, and all the (blue) children of be clockwise from ao become 
children of a^. It is important to note that in order for the peak to valley move 
to be possible, a must be a leaf of the blue tree. Figure [16] shows how the red 
and green trees might look. Using the vertex condition (Figure lai) for b, going 
clockwise around b there are fci > blue edges coming into b after at, followed by 
6's red edge, then there are k2 > green edges coming in, and finally 6 has its blue 
edge be. As above, r„ is the head of v's red edge. We will go through 5 stages in 
the canonical path. 
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Figure 16: Canonical path to move a peak down to a valley in the blue tree. 



Stage 1. If there are green edges coming into b, let d be the tail of the last incoming 
green edge; if b has no green edges coming in, then set d — During the first 
stage of the canonical path, &'s red edge must move from rt, to d (if 7^ c?) by 
making red/green swaps. 

Stage 2. Let / be the last blue incoming edge to b. Now we can move fci + 1 blue 
edges pointing to b (namely, all blue edges between and /I, inclusive) to point 
to c. Do this sequentially and counterclockwise, using &'s red edge to swap. See 
Figure 16 |b-c). 

Stage 3. Vertex a may have green edges entering between the edges or^ and a^. 
Let g be the tail of a's last incoming green edge, if it exists; otherwise, g — Va- If 
9 '^aj we move a's red edge from to g. 

Stage 4. Using a's red edge, we can move the fci blue edges from c back to a, 
correcting the blue tree. 

Stage 5. To repair the red tree, first we need to move a's red edge to point to 
(Stage 5a). Notice that Vs red edge now points toward a. Thus we need to pass 
a red edge from a to (Stage 5b). Since a is a leaf of the blue tree, this can be 
done using a sequence of red/green swaps. 

Given a transition (Z, W) oi M.ef we must upper bound the number of canon- 
ical paths 7x,y that may use this edge. If {Z, W) is in stage 1, then it moves &'s 
red edge, so we know 6. To get back to X, move 6's red edge counter-clockwise 
until it can't make any more red/green swaps. Hence we don't need to remember 
any vertices. Given X and 6, we can find Y since we know which peak to flip down. 
If {Z, W) is in stage 2, then it moves some blue edge from b to c, so we know both 
b and c. We must record /, but then we can get back to the last configuration in 
stage 1, and using b we can get back to X. If {Z, W) is in stage 3, then it moves 
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a's red edge, so we know a. We need to remember to get back to stage 2, and 
remember / to get back to stage 1. Given a we know h, so we can then get back 
to X. If (Z, W) is in stage 4, then it moves a's red edge, so we know a. To get 
back to stage 3, just move a's red edge as far counter-clockwise as possible. Hence 
to get back to X we need to record Tq and /. If (Z, W) is in stage 5a, then we are 
moving a's red edge, so we know a. Rotate it as far counter-clockwise as possible 
to find riy. To recover y, we just need to know so we can pass an edge from a 
to Ta- If (Z, W) is in stage 5b, then we need to remember a and r^. Move a's red 
edge all the way to to obtain Y . In each of the five stages we need to record a 
maximum of 2 vertices. This implies that in this case there are O(n^) canonical 
paths which use any edge {Z, W). 

We have shown that in each of the four cases above there is a maximum of 
0{n^) canonical paths which use any edge {Z,W). If the move of A4dk affects 
both the top and bottom paths, we can think of this move as two moves, each of 
which affects only the top or bottom path; hence, we concatenate the paths for 
each of those moves. Therefore, if we record a bit to decide if the move of AioK 
affects both the top and bottom paths, as well as a bit to decide which of the stages 
we are in, this implies that across all cases there is a maximum of O(n^) canonical 
paths which use any edge {Z,W). Notice that the maximum length of any path 
jXY is 0{n^). We can now upper bound the quantity A which is needed to apply 
Theorem [l2l as follows: 

A ^ max I ——4 7 V hxYW(X)P'(X,Y) 

< max i 4n V i < 0(n«). 

{z,w}eE{P} ^ 2n \ 

Moreover, we can bound tt* as follows: 



1 1 

xFi C„+2C„ - - 30^ 



TT* = min ■k{X) — — — —7, — > 



Applying Theorems 11 and 12 , we get the following 




n log 30 — log e g 



-log(2e) 
0{n^°\og{n/e)) 



n \og{n/e) 



Therefore M.ef is an efficient sampling algorithm for sampling from the set of all 
3-orientations over any triangulation on n internal vertices. 
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5 Concluding Remarks 



Several questions remain open. The complexity of enumerating Eulerian orienta- 
tions in planar graphs of bounded degree is one of the foremost, as raised by [ 15] . 
Extending our fast mixing result to triangulations with larger degrees is a natural 
open problem; perhaps there is an alternate local chain which can sample efficiently 
from the set of 3-orientations corresponding to any fixed triangulation, without re- 
course to the bipartite perfect matching sampler of |18j . Finally, our slow mixing 
example involves vertices of degree ^l{n), and it would be of interest to find other 
constructions with constant maximum degree. 

Acknowledgments. The last author thanks Stefan Felsner for introducing him 
to the problem of sampling 3-orientations and providing useful links to literature. 
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